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Abstract. Dielectrophoresis (DEP), the motion of polarizable particles in non-uniform electric fields, has become an 
important tool for the transport, separation, and characterization of microparticles in biomedical and nanoelectronics 
research. In this article we present, to our knowledge, the first molecular dynamics simulations of DEP of nanometer- 
sized colloidal particles. We introduce a simplified model for polarizable nanoparticles, consisting of a large charged 
macroion and oppositely charged microions, in an explicit solvent. The model is then used to study DEP motion of the 
particle at different combinations of temperature and electric field strength. In accord with linear response theory, the 
particle drift velocities are shown to be proportional to the DEP force. Analysis of the colloid DEP mobility shows 
a clear time dependence, demonstrating the variation of friction under non-equilibrium. The time dependence of the 
mobility further results in an apparent weak variation of the DEP displacements with temperature. 



PACS. 82.20.Wt Computational modeling; simulation - 61.25.Hq Macromolecular and polymer solutions; polymer 
melts; swelling - 82.70.Dd Colloids 



1 Introduction 

' A polarizable particle in a non-uniform electric field is set 
in motion due to the coupling between the electric field 
gradient and the induced polarization. This effect was named 

' dielectrophoresis (DEP) by Pohl in his pioneering studies in 
the 1950's flT\. DEP can be understood by considering the 
simple picture of a dipole in a non-uniform electric field. The 
dipole will be at least partly oriented in the direction of the 
field gradient and hence one end of the dipole will experience 
a stronger electric field strength than the other. This results in a 
non-zero net force and the dipole is set in motion. 

If the polarizable particle is smaller than the characteristic 
length at which the electric field E changes appreciably, the 
DEP force affecting the particle can be approximated as a 
coupling between its induced dipole moment p and the field, 

-Fdep - (p • V)£;. (1) 

For an ideal dielectric the induced dipole moment is linearly 
proportional to the electric field, p = apE, where ap is the 
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total effective polarizability of the particle that depends on the 
dielectric properties of the particle and the suspending medium. 
Combining this relation with Eq. one obtains a general 
equation for the DEP force, 

-Fdep = ^ap\/E^. (2) 

Two important observations can be made from Eq. (|2j- 
First, the direction of the DEP force is parallel to the gradient 
of the electric field squared, but does not depend on the actual 
polarity of the field. Hence, the same DEP effect can be 
achieved with both DC and AC electric fields. Second, the 
direction of the force also depends on the polarizability of the 
particle. If the particle is more polarizable than the medium, it 
will be attracted to regions of higher field strength. If, on the 
other hand, the particle is less polarizable, it will be repelled 
from the high-field regions. The first case is commonly referred 
to as positive DEP and the latter as negative DEP. 

With modern microfabrication techniques based on photo- 
and electronlithographies it is possible to routinely manu- 
facture smooth electrode systems in the sub-micron scale. 
Such miniaturized structures can generate high electric field 
strengths and field gradients with low applied voltages, thus 
minimizing the unwanted effects of electrothermally induced 
motion of the suspending fluid Microelectrode structures 
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employing DEP have been used, e.g., in the self-assembly 
of micrometer-sized silicon resistors |4| and trapping of 
conducting nanoparticles |5 6|, as well as alignment and 
assembly of metallic nanowires L7.8J and carbon nanotubes 1 9 
1101 . The precision obtained in these studies and the fact that 
many particles can be positioned in the electrode systems in 
a parallel fashion show that DEP is a very promising tool for 
the assembly of nanoelectronics devices on a larger scale (for 
recent reviews on the use of DEP in nanotechnology, see Refs. 
1 11 12 1). Since DEP is a non-invasive and non-destructive 
technique, it has also been used to manipulate biologically 
relevant microparticles in various applications L13. 14.15. 16t 
I17I18I19 20|. 

In comparison to standard electrophoresis (i.e., particle 
motion due to the coupling between an applied electric field 
and a non-zero net charge), DEP has the disadvantage that 
the polarization forces involved are typically quite weak. 
In general, efficient particle manipulation in microelectrode 
systems requires balancing DEP against other forces present, 
such as viscous, buoyancy, and electrohydrodynamic forces 
i3i . Clever experimental setups have indeed been devised to 
trap particles of the order of only a few nm in diameter 
fT(3|. But while many of the hindrances can be overcome 
with appropriate choices for suspending media and electric 
field configurations, a fundamental limit on the precision and 
efficiency of DEP manipulation is set by the thermal motion of 
the particles. For the case of stable particle trapping over long 
periods in time, the DEP force clearly has to be much stronger 
than the thermal noise from the particle-solvent interactions. 
To this end, reasonable criteria of the trapping conditions for 
isolated particles with idealized geometries can be formulated 
by using simplified hydrodynamics |14 16). However, in the 
event of particle aggregation the situation becomes much more 
complicated li21i22J . 



Computational modeling has previously been used to study 
diverse electrokinetic phenomena, such as electrophoresis 1 23 
^A\, overcharging and like-charge attraction |25 26 27 1, and 
structure formation in electrorheological fluids |28 29|. It 
is then natural to ask whether information on DEP motion 
of nanoparticles could be obtained from dynamic computer 
simulations employing simplified, but still physically sound 
models for particles and their interactions. The motion of 
nanoscopic particles and particle aggregates in electrode sys- 
tems of realistic dimensions is essentially a multi-scale prob- 
lem, with important processes taking place at times ranging 
from picoseconds (dynamic processes on the molecular scale) 
to microseconds (larger aggregate formation), and all the way 
to seconds (particle concentration in specific parts of the 
device). At present there exists a variety of methods to tackle 
multi-scale problems at different levels of accuracy 13013 11321 
|33|. 

As the initial step to computational modeling of DEP, in 
this article we present, to our knowledge, the first molecular 
dynamics (MD) simulations of dielectrophoresis. The goal of 
the present work is to model the DEP of a single colloidal 
particle at varying combinations of temperature and applied 
electric field strength. We verify that the model exhibits linear 
drift velocity to DEP force relation and further study the 



particle DEP mobility and the temperature dependence of the 
motion. 

This article is organized as follows. In Section |2j we will 
outline in detail and validate the nanoparticle DEP model used 
in the present work, and further examine the extent of finite- 
size effects in our model system. Section|3lthen comprises the 
simulation results on the single-particle DEP transport. Section 
Incloses the article with the concluding remarks. 



2 Simulation method 

2.1 General features of the model 

We have used classical MD simulations to model the dielec- 
trophoresis of a single charge-neutral colloidal particle in a 
non-conducting solvent. As a generic model for a polarizable 
particle we used a system consisting of a large sphere, 
carrying a charge qmac on which A^mic small particles with 
a charge gmic were electrostatically bound (see Sec. \2.2l . In 
the remainder of the text, we adopt the following terminology: 
the large sphere will be called macroion, and the small bound 
particles microions. The macroion-microion complex will be 
called colloidal particle or simply colloid. As the solvent 
particles are smaller than the colloid in linear dimension only 
by an order of magnitude, the colloid can be viewed as a 
nanoparticle immersed in a molecular solvent. 

For a proper description of colloid-solvent interactions 
under non-equilibrium conditions, the suspending medium in 
which the colloid was immersed was simulated with explicit 
neutral particles. Excluded volume interactions between all 
particle types were modeled using a truncated and shifted 12-6 
Lennard- Jones potential (also called the WCA potential) I34l. 
In order to take into account the size of the macroion, an 
additional hard-core radius j q was included in the denominator 
of the repulsive and attractive power terms |35J, 



^wcA = 4e 



(3) 



where rij ~ \rj — rj| is the center-to-center distance between 
the interacting particles i and j (r i being the position vector of 
particle i). The cutoff distance of the potential was tq + 
with ro = 4.5(7 for interaction pairs involving the macroion and 
T'o = otherwise. The characteristic parameters a and e of the 
potential were chosen as the basic units of length and energy, 
respectively. With m further designating the unit of mass, the 
unit of time is then defined as r = a^/m/e. 

The simulation cell was cubic with periodic boundary 
conditions and a side length — 35cr. The solvent particles 
and microions were assigned a mass Im and the solvent 
number density was chosen as ps = 0.3 . For the macroion, 
the mass was chosen to correspond to the mass of bulk solvent 
equal to its volume, giving a value of 155r7i. Suitable estimates 
for effective particle sizes could be obtained from the distances 
of closest approach between different particles at temperature 
T, using the relation <?vvcA = k^T. Our choice of ro then 
gave effective particle radii roughly Scr for the macroions and 
(J 1 2 for the microions and solvent particles. Although particle 
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sizes defined this way depend on the system temperature, the 
steepness of the WCA potential slope ensured that the effective 
particle radii did not vary significantly in the temperature range 
used in this study. 



appreciably. This type of an electric field geometry is a good 
approximation for many experimental DEP setups |12|, and 
was thus a reasonable choice for a model system. The DEP 
force arising from the non-uniformity of the field is 



2.2 Electrostatic interactions 

The electrostatic interactions between the ions were calculated 
directly from Coulomb's law. The energy associated with the 
interaction of particles i and j was then given by 



DEP 



(4) 



where qi is the charge of particle i, eq is the vacuum 
permittivity, and the relative permittivity of the medium. 
To ensure that the microions stay bound to the macroion 
surface we chose the regime of strong Coulomb coupling by 
setting (47reoer)^^ = 56 ecr/e^, independent of temperature 
and constant in the entire simulation cell. We did not employ 
any type of lattice summation or multipole methods as this 
would be problematic due to the non-uniformity of the external 
electric field (see Sec. l2.3t . 

The charge ^mac = — 20e, where e is the elementary charge, 
was assigned to the interaction center of the macroion. To 
prevent an electrophoretic component in the external electric 
field-induced motion, overall charge neutrality was imposed 
on the macroion-microion system, gmac + ^mic9mic = 0. The 
number and charge of the microions were then chosen as iV,„ic 
= 10 and qniic = +2e, respectively. 

At a first glance, the choice of spatially uniform value of 
may seem controversial, as the magnitude and direction of the 
DEP force generally result from the mismatch of the complex 
permittivities of the manipulated particles and the suspending 
medium. A purely practical reason for this simplification 
is the fact that taking into account the induced surface 
charges of polarization (image charges) on the macroion 
l36i is computationally extremely intensive, making dynamic 
simulations, such as the ones presented here, dramatically 
more time-consuming. However, considering the picture of a 
single nanoparticle in solution, it is not even meaningful to 
talk about its permittivity. With our generic macroion-microion 
system, the polarization of the colloid originates simply from 
the redistribution of the microions on the macroion surface. 



2.3 External field 

For the modeling of DEP, it is preferable to use an electric 
field geometry which is simple to realize and has a large 
field gradient in order to produce appreciable dielectrophoretic 
motion. In this work we have used a spherically symmetric DC 
electric field. 



E{R) = Eo[^ 



(5) 



where Eg is the electric field strength at the initial radial 
position Rq of the colloid, and is the radial unit vector. 
As the colloid moves in the field, Rq further determines the 
characteristic length over which the field strength changes 



El 

Uq 



Rq 



(6) 



as results from Eq. Since our model colloid is always more 
polarizable than the surrounding medium, the simulations cor- 
responded to the case of positive DEP. For all the simulations 
presented in this article we used Rq = 1500ct. 

The electric field strength affecting a charged particle was 
calculated by using an additional coordinate system, with 
the origin set at the force center of the electric field. The 
vector R was used to denote the position of the center of the 
microscopic simulation cell in the electric field coordinates. 
The position of a charged particle i was then r[ = R + Ti. 
(In the following, all coordinates marked with prime will refer 
to the electric field frame of reference.) We emphasize that the 
only effect of assigning the additional electric field coordinates 
for the charged particles was to take into account the spatial 
dependence of the electric field and its gradient. During the 
course of the simulations the vector R was updated, and hence 
all as well, according to the colloid displacements in the 
simulation cell (see below in Sec. l2.4> . 

In the absence of an external electric field the microions 
are evenly distributed on the surface of the macroion due 
to their mutual electrostatic repulsions. When an external 
field is applied, polarization of the system is obtained with 
the redistribution of the microions. In order to calculate the 
DEP force affecting the colloidal particle, its polarizability ap 
needed to be determined. This was done by calculating the 
induced dipole moment of the particle, p = g^r^, where the 
summation is carried out over the ions, in the following way. 
The colloid was put under a selected electric field strength and 
the microions were allowed to freely redistribute themselves 
on the macroion surface. The total kinetic energy of the system 
was very slowly quenched to zero, resulting in a minimum 
energy state of the polarized colloid. Note that this type of 
MD energy minimization cannot rigorously guarantee that the 
obtained ion configuration is the true global energy minimum 
state of the system |27|. However, several independent energy 
minimization simulations showed only negligible differences 
in the total energy and thus the obtained values were assumed 
to he very close to the true energy minimum. 

As expected, the dipole moment of the macroion-microion 
system was observed to be linearly proportional to the electric 
field strength in a wide range of Eq. Deviations from the 
linearity occured at very low field strengths due to microion- 
microion electrostatic correlations and at extremely high field 
strengths (above 17 e/ecr, which is nearly half the electric 
field strength due to the macroion bare charge at a distance 
of ro + cr), as the microions were stripped from the macroion 
surface. A linear fit to the data gave the colloid a total 
polarizability ap = 4.94 ±0.16 e^a^/e. Comparing the DEP 
forces calculated with fit value of ap and Eq. (|2j to the 
ones obtained directly from the simulations showed excellent 
agreement, thus vahdating the dipole approximation. 
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2.4 DEP simulations 

The DEP simulations were carried out as follows. The polar- 
ized colloid, in its minimum energy state, was set at the center 
of the simulation cell and the remaining free volume was filled 
with randomly distributed solvent particles. In order to avoid 
strong repulsion forces due to excluded volume overlaps, it 
was ensured that all interparticle distances were larger than 
the WCA potential cutoff distance The system 

was then thermalized for 85t at a selected temperature (T = 
0.03 — 2.4 e/fce) by coupling all the particles to the thermostat 
by Berendsen et al. \ 37 1 with a coupling time constant of 0. 14t. 

The range of temperatures used in this study may sound 
broad in view of experiments. However, just like the electric 
field strength in our model is a measure of the DEP force, the 
temperature, in turn, is simply a measure of the random thermal 
forces and friction affecting the colloid. Note that since there 
were no attractive forces between the solvent particles and, 
on the other hand, the solvent density was much lower than 
required for hard-sphere crystallization of entropic origin, the 
solvent remained in a fluid state even at the lowest simulation 
temperatures. 

After thermalization the thermostat coupling was retained 
only within a distance of la from the cell borders for the 
rest of the simulation. The border scaling was maintained in 
order to dissipate the heat produced from the work done by the 
electric field on the colloid, without affecting its DEP motion. 
A somewhat similar thermostat coupling scheme was also used 
by Tanaka and Grosberg in their simulations of nanoparticle 
electrophoresis 1231 . In those studies it was reported that the 
thermostat at the cell borders serves only a minor purpose, 
as the amount of Joule heat produced during the simulations 
was very small. Our test simulations without the cell border 
thermostat also showed only negligible increases in the tem- 
perature. We nevertheless retained the border thermostat in the 
DEP simulations to dissipate any excess heat in the system. 

The equations of motion were integrated using the velocity 
Verlet algorithm 1 38 1. For proper energy conservation, integra- 
tion time steps between 0.002t and 0.004t, depending on the 
system temperature, were used. The total simulation length, 
including the system thermalization time, was 5785r. During 
the simulations the colloid, driven by the DEP force, could 
move distances corresponding up to a few simulation cell side 
lengths (see below in Sec. l2.5> . In order to avoid artifacts due 
to the direct calculation of the electrostatic interactions and the 
colloid entering the thermal scaling region at the simulation 
cell borders, after each time step the origin of the simulation 
cell was set at the interaction center of the macroion. In 
practice, this means that the macroion itself was not displaced 
in the simulation cell, but the positions of all the microions 
and solvent particles were shifted instead. The position of the 
center of the simulation cell in the electric field coordinates, 
R, was also shifted according to the macroion motion, as the 
electric field strength and its gradient changed with the colloid 
displacements. 

For each combination of electric field strength and temper- 
ature, a total of 20 - 80 simulations were carried out in order 
to obtain good statistics. Larger numbers of simulations were 
required at higher temperatures due to thermal noise. 



Although the simulations were carried out in the regime of 
strong Coulomb coupling (cf. Sec. 12. 2K it is still valid to ask 
whether the microions could be stripped from the macroion 
surface by random collisions with the solvent particles. In 
that case the colloid mobility would be enhanced due to an 
electrophoretic component, as the particle would possess a 
non-zero net charge. We determined the minimum microion 
binding energy for the case of the highest field strength, Eq 
= 2.70 e/ecr, while the colloidal particle was held fixed at its 
initial position in the electric field. This analysis showed that 
even at the highest temperature used in the simulations, T = 2.4 
e/fce, the ratio of the minimum microion electrostatic binding 
energy to the thermal energy fceT was roughly 24. Hence, the 
macroion-microion system could be considered stable in view 
of collisions with the solvent particles. All the DEP simulations 
were further checked for possible microion detachments, but no 
such events were found. 

If the electric field-induced colloid motion in the simula- 
tions was truly due to DEP, similar deterministic motion should 
be obtained regardless of the actual polarity of the electric 
field. However, if the direction of the colloid motion would be 
reversed by reversing the field polarity, the motion would be 
in fact due to electrophoresis. To make the final validation that 
the model presented here really produces DEP motion of the 
colloidal particle, 10 independent test simulations were carried 
out with both positive and negative field polarities and |i?o| = 
0.96 e/ea. In order to facilitate the comparison, the solvent 
thermal noise affecting the colloid was minimized by selecting 
a very low temperature, T = 0.03 e/fce. The simulations were 
otherwise identical in the initial conditions, except that for 
the case of negative field polarity, the microion configurations 
on the macroion were rotated half a turn with respect to the 
radial direction. Thus, for both field polarities the colloid was 
close to their minimum energy state, and delayed coupling of 
the colloid to the electric field due to microion distribution 
relaxations did not hinder comparison between the two series 
of simulations. 

In all the cases, regardless of the field polarity, the colloid 
moved in the direction of WE'^, as was expected for positive 
DEP. The average colloid radial displacements for the positive 
and negative field polarities were (— 14.0±2.5)<t and (— 12.5± 
2.5)cr, respectively. Considering the number of test simulations 
carried out, it can be stated that the agreement is very good. 

2.5 Finite-size effects 

Before proceeding to the actual DEP simulations, the extent 
of finite-size effects in our model should be discussed. Since 
in general hydrodynamic effects, mediated by the solvent, are 
long-ranged f24"39|, regardless of the simulation cell size 
used there always remains some contribution from interactions 
between the colloid and its periodic images. This imposes 
restrictions on the minimum size of the solvent shell used in the 
modeling that in principle need to be resolved for each specific 
system. 

Probably the best point of comparison in the literature for 
our model system is the one used by Tanaka and Grosberg 
in their MD simulations of nanoparticle electrophoresis |23|. 
The authors used a similar WCA solvent and also the same 
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Fig. 1. Colloid average radial displacements as a function of the 
simulation cell side length Lq. The error bars in the figure correspond 
to the standard deviations of the radial displacements. 



solvent number density as in our model. The colloid radius in 
their simulations varied between 3(T and la. From tests on the 
simulation cell size dependence the authors concluded that a 
value of Lq = 2Qa was sufficient to yield statistically reliable 
electrophoretic drift velocities. The actual electrophoresis 
simulations were then carried out with Lq = 32cr. An important 
difference between our simulation system and the one used 
by Tanaka and Grosberg is that the latter included co- and 
counterions in the solvent, which the authors stated to shield 
hydrodynamic interactions. It should be also noted that in prin- 
ciple the thermostat coupling at the simulation cell borders we 
have used could affect the colloid hydrodynamic interactions. 
However, as the thermostat coupling was quite weak, this effect 
should be very small or even negligible. This is corroborated by 
the tests of Tanaka and Grosberg on the cell border thermostat 
and nanoparticle electrophoretic drift velocities |23 1. 

We tested several simulation cell sizes in order to find a 
value of Lq small enough for good computational efficiency, 
but still large enough to result in reasonably small finite-size 
effects. As the measure for an adequate value of Lq we used 
the colloid average radial displacements, (AR), in reference 
to the electric field coordinates, at Eq = 0.96 and 1.91 e/ecr. 
Like in the tests of opposite electric field polarities (Sec. l2.4L 
we used a temperature of 0.03 e/Zcb to minimize statistical 
fluctuations. Figure shows the results from several tests 
employing simulation cell side lengths between 20a and 42.5(7 
with 2239 and 22728 solvent particles, respectively. 

For both Eq used the absolute values of the colloid average 
radial displacements (AR) increased with decreasing Lq. This 
is in contrast with the well-known enhancement of hydrody- 
namic friction, proportional to L^^ II24I39II . However, standard 
hydrodynamic considerations assume an incompressible fluid, 
which is not the case with our system. In the simulations of 
nanoparticle electrophoresis, Tanaka and Grosberg also ob- 
served 1 23 1 increased particle drift velocities for smaller values 
of Lq. Although in their case the solvent size dependence of the 
particle velocity was not smooth as observed here, the effect 
could be due to the general properties of the low-density WCA 
solvent. 

A measure of the relative efficiencies of momentum and 
mass transport in a fluid at equihbrium is given by the Schmidt 



Sc = 



PmDs 



(7) 



where 77, p„i and Dg are the fluid shear viscosity, mass density, 
and self-diffusion coefficient, respectively. In real liquids Sc is 
typically of the order of lO'^ and thus momentum propagates 
much faster than mass in the medium. However, for our 
model solvent, simple estimates based on the theory of hard 
sphere fluids |40| result in values of 5c w 3 - 5 (depending 
on the temperature). This suggests that the apparent friction 
on the colloid under non-equilibrium conditions is strongly 
connected to the number of solvent particles and hence, Lq. 
Considering the fact that momentum is not conserved in our 
DEP simulations (due to the work done by the electric field 
on the colloidal particle and also, to some extent, due to the 
temperature scaling only at the simulation cell borders), the 
exact coupling between the solvent flow field and the colloid 
dynamics becomes hard to assess with analytical means. 

Based on the tests, we selected a value of Lq — 'iba, with 
12687 solvent particles, for the actual DEP simulations. This 
value of Lq still seems to produce larger radial displacements 
than would be obtained for a system with a much higher 
number of solvent particles. However, since all the simulations 
were carried out with the same system size, and our purpose is 
to provide qualitative rather than quantitative insight on DEP, 
we do not expect finite-size effects to affect the conclusions of 
the present work. 



3 Results 

In this Section we present the analysis of the DEP simulations. 
In all the cases the initial position of the colloidal particle 
was set at a distance _Ro = 1500(t on the x'-axis. Hence, at 
the start of the simulations the DEP force drove the colloid in 
the negative a;' -direction, while motion in the y'z'-plane was 
perpendicular to the DEP force. The DEP force affecting the 
colloid was varied using four different electric field strengths 
with = 0.96, 1.35, 1.91, and 2.70 e/ecr. Since in our electric 
field geometry Fdepoc E"^, the ratios of the initial DEP forces 
were approximately 1:2:4:8. The simulations were carried out 
at temperatures between 0.03 and 2.4 e/fce- The margins of 
eiTor for the quantities presented below are estimated from the 
error of the mean. 

In view of DEP manipulation in the nanoscale, it is 
instructive to compare the DEP force to the thermal noise of 
the suiTounding solvent. For the latter, a characteristic force is 
given by the ratio of the thermal energy to the linear size of the 
particle in question. 



thermal 



2ao ' 



(8) 



where in our case ao ~ ro + (t/2 is the radius of the spherical 
macroion. A dimensionless parameter describing the ratio of 
the thermal noise to the DEP force, is then defined as 







DEP 



thermal 



knT 



a^anWE^l' 



(9) 
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Fig. 2. Average colloid radial displacements, as a function of Eq, 
at the end of the simulations for three different temperatures. The 
solid lines are the results of the regression analysis between the radial 
displacements and the spatially variant DEP force (see text). 



At small values of Obep the DEP force dominates over the 
thermal motion of the colloid, whereas at large values the 
colloid motion becomes increasingly random. In the range of 
electric field strengths and temperatures used in this study the 
values of Ouep varied nearly three orders of magnitude, from 
0.07 to 40. However, there is an important point to note here. 
As the DEP force is strongly dependent on the radial position 
in the electric field (Fdepoc R^^), the value of Ouep is not 
fixed. For example, a colloid that has traversed a distance of 
lOOcr radially inward in the electric field experiences a DEP 
force that is 41% larger than the force at its initial position, 
Rq = 1500(7. This further entails that at sufficiently long times 
particle distributions in the radial direction become asymmetric 
with a negative value of skewness. 

It can be questioned whether the DEP force originating 
from the polarization of the macroion-microion complex is 
the same at all temperatures. As the solvent particles collided 
with the microions, the dipole moment of the polarized colloid 
fluctuated around the value corresponding to the minimum 
energy state. However, even at the highest temperatures used, 
these deviations were not strong enough to divert the mean 
orientation of the dipole moment from the one induced by 
the electric field. For all cases the characteristic fluctuation 
time of the dipole moment was very short in comparison with 
the length of the simulation. It was then plausible to take 
the average force affecting the colloid as the one given by 
Eq. (|6jl. Test simulations with a spatially dependent, but non- 
fluctuating, artificial force affecting the colloid showed similar 
radial displacements as in the DEP case (data not shown), 
further coiToborating the approximation. 



3.1 DEP at constant temperature 

The main results of our simulations were the colloid radial 
displacements in the electric field, i.e., motion in the direction 
parallel to the DEP force. We first consider series of simula- 
tions where the system temperature was kept constant and Eq 
was varied. The average radial displacements of the colloid, 
(AR), for three different temperatures are shown in Fig.|2| The 



10' 

5 

I 10 

5 

2 
1 



♦ 

A 
□ 



4 

■A 



£„ = 2.70 elexT 
£„ = 1 .91 
Eij = 1 .35 
Eij = 0.96 



' ' "I 

£]■■■ 



♦ 



...[] 



10" 



T [elks) 



Fig. 3. Absolute values of the average colloid radial displacements 
for the different values of Eo as a function of temperature. 



increase of the radial displacements with the increase of Eq and 
the decrease of T is evident. 

The ratios of the DEP drift velocities, (udep) = {AR)/At, 
and the colloid thermal velocities at the respective temperatures 
were of the order of 0.01 - 1 . At these low values the DEP drift 
velocities should be linearly proportional to the DEP force. The 
average radial displacement over a time interval At is then 



(AR) 



AL 



(10) 



where ^ is the friction factor of the solvent. Regression analysis 
between the particle radial displacements and the DEP force, 
(AR) cx [i^DEp]'^, justified this assumption. However, for an 
accurate analysis it was required to taken into account the fact 
that the particle in our simulation is under spatially variant 
force field: As the radial position of the colloid changes, so 
does the magnitude of the electric field and its gradient. Hence, 
the radial displacements were compared to the average DEP 
force experienced by the colloid, as calculated from the particle 
trajectories. The resulting values of v were, within statistical 
eiTors, equal to unity (see Fig.|2j. 

On the other hand, simply taking the value of the DEP force 
at the start of the simulation, i^DEp(^o), resulted in values of 
1/ slightly higher than one. 



3.2 Temperature dependence and DEP mobility 

The data in Fig. |2] indicate that the dependence of the DEP 
displacements on temperature is quite weak. We now analyze 
this observation in more detail. 

Figure |3] shows the colloid radial displacements, as a 
function of temperature, for four series of simulations corre- 
sponding to different values of Eq. Note that the data presented 
are the absolute values of the radial displacements (in order to 
make the full logarithmic plot) but all the real values of the 
displacements were in fact negative, i.e., in the direction of the 
DEP force. The data show a clear power law behavior, although 
for the two lower values of Eq, at the highest temperatures, the 
diffusive motion of the colloid contributes significantly to the 
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scatter in the DEP displacements. This is characterized by high 
values of the parameter Ouep ~ 8 — 40. 

As the linear response was shown to hold above, a 
straightforward assumption is that the friction factor ^ should 
be linearly proportional to the solvent shear viscosity ?/. 
An estimate for the temperature dependence of our solvent 
viscosity can be obtained from the well-known Enskog ex- 
pression for hard-sphere fluids |!40|. It has been shown to give 
good agreement with MD calculations employing low-density 
WCA potential solvents |41| (such as the one in our case) 
and further even qualitatively describe much more complex 
systems where one would not necessarily expect the hard- 
sphere approximation to hold 1421 . For the Enskog viscosity 
the temperature dependence is 77 c>c VT. Although one could 
assume even a qualitative agreement with our low-density 
WCA solvent friction and the one given by the Enskog theory, 
this is clearly not the case (cf. Fig.|3}- What is then the reason 
for the scaling of the DEP displacements shown in Fig.lsj? 

Analogous to the case of electrophoresis under uniform E, 
it is possible to define DEP mobility 1 17 1, /^dep, that yields the 
particle drift velocity due to the electric field, 

(^^DEP) - MDEPV£:2. (J J) 

Combining Eqs. (|2j, dlOt . and il 1> yields the general form 



a 



MDEP = 777- 



p_ 

2i- 



(12) 



However, care should be taken here. As mentioned above, 
contrary to the case where true steady-state condition can be 
attained, the value of VE^ changes constantly with particle 
motion. Calculating /idep from large particle displacements 
and simply using the initial value of VE^ can lead to significant 
errors. 

Using Eq. (II 1> we calculated the colloid mobilities from 
the particle trajectories and the average values of VE'^, over a 
given time interval At from the beginning of the simulations. 
The results of these calculations for two different temperatures 
are shown in Figs.|3and|5l For clarity, we only show the results 
for two different values of i?o in each plot. The data for all Eq, 
at each respective temperature, were consistent with each other 
Also note that for lower values of Eq the data are somewhat 
noisier due to the stronger diffusive motion. 

The main observation here is that the values of /idep, 
and hence the values of the friction factor ^ (cf. Eq. (II \V ). 
vary with time. Thus, there is a clear coupling of the non- 
equilibrium motion of the colloid to the solvent velocity field. 
Although a thorough analysis of this effect is challenging, 
it is not unreasonable to assume that the characteristic time 
scale for the evolution of the non-equilibrium transport (with 
a constantly changing friction factor) would be different at 
different temperatures. Indeed, a further study of the apparent 
temperature dependence of the DEP displacements (such as 
shown in Fig. |3l for the final radial displacements in the 
simulations) revealed that while power law expressions could 
be fitted to the displacement data at different points in time, 
the resulting exponents were not constant. We do not want 
to put too much emphasis on the actual values of these data, 
as we were unable to find a satisfactory theoretical way of 
interpreting them. 
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Fig. 4. Colloid mobility /idep as a function of time for Eo = 0.96 and 
1.91 £/eaatT = 0.03 e/fce. 
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Fig. 5. Same as Fig.lHbut for Eq = 1.35 and 2.70 e/ea at T = 1.0 



3.3 Lateral displacements 

Another interesting aspect of the spatial dependence of the 
colloid motion is seen by studying it in the plane perpendicular 
to the electric field gradient at the initial particle position, 
i.e., in the y'z'-plane. At the beginning of the simulation the 
colloidal particle is directed toward the negative x'-direction 
by the DEP force. Deviations from this direction then affect 
the precision of the DEP manipulation, i.e., how accurately the 
final position of the particles can be predetermined. 

To a very good approximation the average displacements 
of the colloid in the y'- and ^'-directions were zero for 
all the combinations of Eq and T. We then calculated the 
displacements in the j/'z'-plane, /\i?iat. For comparison, we 
also determined the collo id tra cer diffusion coefficient D = 
(3.65 ± 0.07) X 10^2 cTy/^J^ at T = 1.0 e/fce in the absence 
of the external electric field. This was done using the Green- 
Kubo relation between the diffusion coefficient and the particle 
velocity autocoiTelation function 1401 . 



(13) 



A total of 10 independent simulation runs were carried out 
with a length of 15800t each. The integrals of the velocity 
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Fig. 6. Lateral {y' z') displacements ARi^t as a function of time at T 
= 1.0 e/fce- The solid line is the prediction obtained from Z\7?iat(t) = 



autocorrelation functions were calculated over time periods 
of 158t with a sampling interval of 0.002t to ensure proper 
convergence. 

A comparison between the average lateral displacements 
for the two extreme cases _Eo = 0.96 and 2.70 e/ea, and 
the root-mean-square displacement predicted by the relation 
^R\a,t= V 4I?< is shown in Fig. |6l In the case of Eq = 
0.96 e/ecr, for which (AR) was small, the colloid average 
lateral displacement is in a good agreement with the zero- 
field diffusion prediction. However, at Eq = 2.70 e/ea a clear 
difference is seen. This indicates that with decreasing radial 
position in the electric field, the high electric field gradient 
produces a steering effect on the colloid. Deviations from 
motion along the x'-axis result in DEP force components also 
in the y'- and z'-directions, and the assumption of independent 
motion in the a;'-direction and, on the other hand, in the y'z'- 
plane breaks down. It is also possible that the coupling of 
the particle motion to the solvent velocity field suppresses the 
lateral motion up to some degree. 

A general relation between the solvent friction and diffu- 
sion coefficient of a particle is given by the Einstein relation 

D^'-f, (14) 

regardless of the actual mechanism causing the friction. From 
the simulated value of D we calculated the friction factor 
^ = 27.4 ± 0.5 y/em/a, coiTesponding to the case of zero 
external electric field. It is interesting to compare this value 
to the time-dependent ^ in the DEP simulations. Inserting the 
value of ^ above into Eq. J12l i results in /idep = (0.090 ± 
0.003) e^cr^/re^ in equilibrium. This is comparable to the 
values of /xdep in the non-equilibrium simulations only at 
very short times, see Fig. |5] Although not conclusive, the 
analysis here implies that at the beginning of the simulation the 
DEP mobility is very close to the value obtained by using an 
equilibrium friction factor. The coupling to the solvent velocity 
field then gradually lowers the solvent friction factor in time 
and the particle motion loses its equilibrium nature. 



4 Concluding Remarks 



We have carried out extensive MD simulations of a spherical 
nanoparticle DEP in a non-conducting solvent. This study is, 
to our knowledge, the first MD study of DEP and serves as a 
benchmark for future work employing more complex colloid, 
solvent, and electric field models. 

The main results of our simulations were the radial 
displacements of the colloidal particle due to the DEP coupling. 
The analysis of these displacements first showed that the 
assumption of linear response was justified for simulation 
series carried out at constant temperatures. That is, the DEP 
drift velocities were Unearly proportional to the DEP force. For 
an accurate analysis it was required to take into account the 
spatial dependence of the electric field, even at the short time 
scale of the present simulations. 

The comparison of simulation series at different temper- 
atures showed a weak apparent temperature dependence. We 
emphasize here again that the magnitude of this dependence 
was observed to vary with the time scale over which the particle 
motion was studied. Furthermore, the particle DEP mobilities, 
defined by Eq. were observed to increase in time. This 
is clearly a manifestation of the non-equilibrium nature of the 
particle motion. The intricate coupling of the particle drifting to 
the solvent velocity field decreased the friction factor (which 
is inversely proportional to the mobility). A rigorous analytical 
treatment of this effect is indeed challenging, and is left as a 
subject for future work. 

Finally, the analysis of the colloid lateral displacements 
showed that while the DEP mobility itself was shown to be 
time-dependent, for the lowest value of Eq the colloid motion 
perpendicular to the initial direction of the DEP force could still 
be well predicted by an equilibrium diffusion coefficient. With 
increased DEP coupling, this no longer applied and the lateral 
motion became coupled to the DEP drifting. Intuitively, this is 
not surprising. However, it is an important effect in determining 
the precision of particle transport, and we are not aware of 
any direct quantification of such steering effects neither in 
experiments nor computational studies of DEP. 

Although thermal motion of the nanoparticles sets lim- 
itations on the DEP forces required for their deterministic 
motion, other hampering effects such as the fluid streaming 
induced by local heating j3| can be suppressed with proper 
experimental set-ups. For these conditions simple dynamic 
particle simulations, such as the ones presented here, can yield 
valuable insight on the conditions, precision, and time scales 
of efficient transport and stable trapping. This is especially 
important with particle sizes decreasing to the nanometer 
range, where the direct observation of the manipulated particles 
and their aggregates becomes very difficult in experiments. 
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